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cases  were  run  and  initiation  behavior  was  plotted  both  as  a function  of  the 
total  energy  of  the  bursting  sphere  and  the  energy  density  of  the  bursting  sphere 
In  addition,  a few  runs  were  made  for  the  ramp  edition  of  energy  to  compare  to 
the  bursting  sphere  results.  The  study  shows  that  as  the  energy  density  of  the 
sphere  ges  becomes  lower  the  total  energy  required  for  initiation  increases 
sharply.  This  is  in  agreement  with  experimental  studies  using  spark  initiation. 
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I.  INTRODUCTION 


Since  the  early  1940' s,  there  liave  been  a number  of  models  for  the 
process  called  detonation.  These  range  from  the  early  one-dimensional 
model  of  Zel'dovich  and  Kompaneets^  to  models  taking  account  of  three-di- 
mensional structure  including  those  of  Barthel  and  Strehlow^  and  Barthel^. 

The  one-dimensional  model  is  well  developed  and  the  three-dimensional 
model,  though  less  well  understood,  yet  provides  much  qualitative  insight 
into  the  mechanism  responsible  for  the  self-sustaining  nature  of  the 
detonation. 

Within  the  scope  of  the  .subject  of  detonation,  one  can  specifically 
focus  attention  on  the  initiation  process.  Initiation  is  usually  considered 
to  include  the  following  event.s.  first,  a shock  wave  travels  through  a 
combustible  mixture,  raising  its  temperature  and  pres.sure.  Next  an  induction 
sone  develops  in  which  chomicul  reaction  begins  to  take  place.  The  temp- 
erature and  pressure  in  this  r.one  lomain  fairly  constant  as  any  heat  release 
is  ai)proximatcly  cancelled  by  a corro-ponding  endothermic  dissociation. 

Tlie  induction  zone  is  characterized  by  the  induction  delay,  also  called  the 
delay  time  to  ignition,  which  Is  simply  the  time  a particle  of  the  mixture 
spends  in  the  induction  zone.  I•■inally,  the  recombination  or  heat  release 
region  develops.  This  of  cour^o  yields  the  observed  e.\othermic  nature  of 
a detonation.  As  can  be  seen  this  on>:ire  process  is  non-steady. 

It  would  be  instructive  to  obtain  a model  for  titis  initiation  process 
witich  will  yield  the  qualitative  m'havior  which  re.'tults  as  one  varies 
input  parameters  such  as  tomperatiu'e  and  shock  overpressure.  This  then  is 
tlte  scope  of  this  paper.  A model  for  the  initiation  of  shock  induced  deto- 
nation is  developovl  for  a spherical  flow  field.  Two  different  types  of 
nou-iJe.il  flow  fields  are  tlten  gen»ruted  in  whieh  tin  initiation  may  occur. 


These  include  flows  generated  by  massless,  higli  pressure,  bursting 
spherical  vessels  and  £lov;s  resulting  from  finite  ratv  depositions  of 
energy  at  a spherical  region  in  space.  Tor  all  these  flow  fields,  the 
model  then  yields  information  cn  delay  time  to  ignition,  minimum 
energy  requirements,  sound  speed  effects  and  explosion  loci. 


II.  BACKGROUND  FLOW 

The  initiation  process  must  occur  in  the  flow  field  under  consider- 
ation. This  flow  field  is  generated  by  a computer  program,  which  was 
obtained  from  A.  K.  Oppenheim.*’  The  CLOUD  program,  us  it  is  calleil  in 
this  paper,  uses  a finite  difference  scheme  to  evaluate  the  fluid  dynam- 
ic conservation  equations  in  onc-dimensionul  Lagrur.gian  form,  These 
include  the  continuity,  momentum  and  energy  equation.''.  A fourth  expres- 
sion, the  equation  of  state  is  also  evaluated.  The  e<iuation  ef  state  for 
this  investigation  takes  the  form  of  tlio  perfect  gas  law. 

The  CLOUD  program  performs  itt.  calculations  with  non-dimensional 
variables.  This  fact  is  crucial  to  lator  computat :on'<  because  it  allows 
ono  to  substitute  in  desired  reference  parameters  after  a given  computer 
run  has  been  made. 


The  indopondont  variable;:  for  tin*  numerical  computation  are  V = t/ll. 


and  U “ r/r^  where  T'and  r denote  uitu'nsiona;  lime  a;vl  radial  position  and 


and  I'y  are  the  reference  p.irameiei  • used  to  ron-d imonsionalite  the 


calculations. 

The  dependent  variables  of  prune  import. uive  are  the  specific  volume 


and  energy,  which  are  given  b/  V = v/v^  and  H = e/T^  respectively,  whore 


again  the  v and  e iutve  dlmenstvms-  anc  the  ftiiiscriptcc  •:ymhu!s  arc  tired  to 
non-dimensional i se  the  caletil-it iv»n *.  Other  depi'iuient  variables  winch  tne 
Ui.OUD  progr.im  calculates  an*  velueif),  pressure  atid  ''seosity.  the  pressure 
at  varlou.s  point?:  in  the  flow  fiel.l  : used  to  deter  ;ne  the  time  of  slioek 
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passage  even  though  the  pressure  does  not  enter  into  initiation  calcula- 
tions. The  viscosity  terra  is  used  to  smooth  out  infinite  gradients  at 
the  shock  location,  as  these  are  unacceptable  in  the  finite-difference 
computational  scheme. 

In  order  to  validate  the  data  concerning  the  flow  field  which  are 
obtained  from  the  CLOUD  program,  three  comparison  runs  were  made  for  the 
case  of  a bursting  sphere.  These  include  sphere  bursts  of  32  atm  air  into 
1 atm  ambient  air  and  15.79  atm  helium  into  1 atm  ambient  air.  These 
runs  are  compared  to  corresponding  runs  made  by  Brode.^  Al.so  a high  pres- 
sut'o  sphere  burst  of  200  atm  air  into  1 atm  ambient  air  is  compared  to  a 
corresponding  run  made  by  Huang  and  Cliou.^  These  data  are  shown  graphically 
on  figs.  1-5.  On  these  plots  the  non-dimensional  shock  ovorpros.sure 
is  shovm  us  a function  of  non-dimtmsionai  energy  scaled  radiu.s  or  shock 
location.  is  Ueterjuinod  from  the  expression 


i* j.  = ’ j.'- 

o 
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where  I*  is  the  value  of  the  shock  ovtnpves.Sure  imd  D is  the  ambient 
s ' 0 

reference  pressure.  K is  obt.'uned  from  the  expression 


u; 

“ r 


where  U is  the  loeutiuu  of  the  shock  and  II  is  eaiUd  the  scaling  radius 
s o 

and  i.^  determined  by  the  following  reUiion,  as  discussed  in  Stvehlow  and 

7 

Uicker. 
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Here  1*^  again  Vepresentf  the  atilnent  nUerenee  pressuie  and  hj  !♦.  the 
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amount  of  energy  which  must  be  added  to  a sphere  at  ambient  conditions  to 
raise  its  temperature  and  pressure  to  the  conditions  necessary  for  the 
burst.  This  energy  addition  is.  assumed  to  occur  instantaneously  and  at  con- 
stant volume,  hj  is  found  from  the  relation. 


1*  - P 

4it  2 0 3 
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where  is  the  pressure  in  the  spliere  just  prior  to  burst,  Jj  is  the 
specific  heat  ratio  of  the  gas  in  the  sphere  and  r^^  is  the  initial  radius 
of  the  sphere.  Then  substituting  eqn.  4 into  eqn.  3 and  the  resulting  form 
of  eqn.  5 into  ecjn.  2,  it  is  found  that 
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It  must  also  be  noted  that  in  Pigs.  1-5  the  first  point  is  not  obtained 
from  the  finite-difference  caicuiatiou  but  is  determjned  analytically 


from  the  shock  wuve-rarefaetkm  contact  surface  ret,uitcttont , as  found  in 


Liepmann  and  Roshko,^* 


1- 


ift-l)  (i«/la>l  IPi/ft  -4 
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where  a refers  to  the  awbienr  sound  speed,  a,  is  the  spheie  sound  speed 
just  prior  to  Inir.st  and  As  the  sinn  ific  heat  rat  io  in  the  acihlent  air. 
Ail  other  quantitliss  are  :he  Mce  as  defined  previousiy.  fqn.  fe  is  itera- 
tively solved  for  the  quantity  «duch  l.s  tl  en  sub.stitoted  lure  eqn.  I 

to  obtain  the  initial  point. 


Also  shown  in  l*igs.  l-d  are  dat.i  obtained  fiwaa  another  nu^.eric.ii 


prugritia.  Tha.s  other  program,  iitown  a . the  WUNik  program  was  oi<tained 

■ . 

5 row  H.  M.  iAternberg.  As  can  l*e  •^een  the  CLOUl"  and  atkliv  data  are  in 
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very  good  agreement.  This  wil.’  then  ho  tni.cu  as  evidence  of  the  accuracy 
and  validity  of  data  obtained  from  tlu  CLOUD  program. 


III.  IGNITION  MODLL 

The  ignition  model  analyiunl  in  this  investigation  is  viewed  as  a 
simple  one-dimensional  model  wliich  will  allow  a parameter  study  of  the 
problem.  General  trends  and  behavior  are  .sought  rather  than  specific 
numerical  values.  With  this  in  mind  then,  the  following  procedure  is 
employed. 

As  given  in  Cusey^^  the  l»>cal  delay  time  tc  ignition  for  a hydrogen- 
oxygen  system  is  given  by  the  vormula, 

T;,  6xp(Fa/RT) 


wliere  A is  called  the  pre- exponential  faetoi'i  (0  J is  the  oxygon  concentra- 
tion, in  the  chemical  activation  energy,  U i;'  the  .mivorsal  ga.s  constant 
f*nd  r is  the  temperature.  Ihe  local  eelay  time  will  s^efer  to  a particular 
particle  or  "cell*'  in  the  floit  field.  Thus  at  lach  instant  of  time,  each 
particle  in  the  flow  will  have  an  instantaneous  value  of  local  delay  time, 
Tj.  Ke  can  then  define  the  variahleV^  to  he  the  fraction  of  total  delay 


time  consumed,  where  the  total  delay  < ime  is  e«p  al  to  the  time  spent  in  the 


induction  2one.  V is  then  given  by  “he  folios ing  expression  taken  from 
Strehlow,^^ 


where  t , refers  to  the  time  of  shock  rcissagu. 
s 

ihe  basic  approach  is  to  iniey.rate  etjn.  5»  sslong  a particle  path 


until  its  value  is  ex|uat  to  initty. 


Tms  will  correspond  to  the  locution 


in  space  and  time  when  the  paiticle  leaves  the  induction  zone  and  ignites, 
The  Lagrangian  approach  is  thus  seen  to  be  very  useful  as  one  must  inte- 
grate eqn.  8 while  following  individual  particles  in  the  flow.  Eqn.  8 
must  be  cast  into  a slightly  different  form  for  use  with  the  CLOUD  program. 
ITiis  will  be  discussed  later. 

It  is  now  necessary  to  specify  the  combustible  mixture  and  the  chemi- 
cal kinetics  to  be  studied.  For  simplicity,  all  computer  runs  were  made 
with  the  expanding  .shock  wave  passing  into  ambient  arr  with  = 1.4,  Air 
is  not  a combustible  mixture  but  it  i.^;  assumed  that  tncre  is  enougii  II2  for  a 
reaction.  In  this  analysis  only  particles  or  cells  which  are  outside  of 
the  pressurized  sphere  or  region  of  energy  addition  are  considered  to  be 
combustible,  llten  as  the  flow  progresses,  a contact  surface  will  always 
separate  these  two  regions.  Cells  in.^ide  of  this  contact  surface  are 
not  followed  as  they  are  not  considert'd  to  bo  combustible.  Cells  outside 
of  the  contact  surl'aco  are  followed  a;'  they  are  considered  to  bo  combustible 
and  hence  are  necessary  for  the  delay  time  calculul.ions.  The  inside  colls 
do  however  serve  the  purpose  of  providing  the  necessary  energy  for  tl^e 
establishment  of  the  shock  wave  which  processes  the  combustible  mixture. 

This  is  illustrated  in  F'ig.  4 for  Ihe  case  of  a bursting  sphere,  .Note 
that  the  outward  travelling  shock  wave  begins  at  the  contact  surface  but 
because  of  the  infinite  gradients  involved,  this  urea  must  InUially  be 
smoothed  out  over  the  .‘•pace  of  two  particles,  krewa  as  fairing  ceils. 

I'or  the  actual  reaction,  .1  simple  hydrogen- o.xygen  system  is  considered. 
Schott  and  Kinsey  “ have  shown  that  the  primary  reactions  involved  are: 


“d'  ‘’.l“ 

11  V u,- 

U * llj- 

uii  * Ur 


11,0  or  Oil 


Oil  ♦ 0 
Oil  f il 
-4  11,0  II 


7 


where  reaction  (ii)  is  the  most  hlgiily  endothermic  reaction  in  the  chain 

(reactions  ii , ili  and  ivj . Fcr  this  system  of  reactions,  Schott  and 
12  13 

Kinsey  and  Strehlow  show  p]ots  of  local  delay  time  times  oxygen 
concentration  vs.  1000  over  the  temperature.  I'rom  these  plots  the  quanti- 
ties E^,  the  activation  energy  and  A,  the  pre-exponential  factor  can  be 
determined.  The  plot  from  Strehlow  is  reproduced  in  IMg.  5.  The  above 
quantities  are  found  in  the  following  way. 


We  have  that  that  the  local  dela>  time  is  given  by  eqn.  7,  repeated 


here 
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that  thi.s  equation  is  of  the  form  y = m.\  + b 
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whore 
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Tlie^  si epus  iuul  intercepts  on  rl\i*  ytrohlow  and  Hclintt  unJ  Kinsey  graphs  ]; 

cun  he  Jetemined  .uid  then  the  vnluc*  of  itid  A fnilow.  These  values  are  | 

siiuun  In  Table  i bi?low. 


Investigator 
Strehlow 
Schott  § Kinsey 


Rto  10 
3541.18 


,?,nA 

’ ' &n  10 

- 10.165 

- 10.647 


(kcal/mole) 

, /sec-moles^ 
liter  ) 

16.19 

6.8418x10"^^ 

18.14 

2.2542x10'^^ 

TABLE  1 


At  this  point  an  arbitrar>  choice  was  made  to  use  the  values  obtained 
from  Strehlow  for  the  initiation  and  delay  time  calculations.  It  was  felt 
that  this  is  justified  by  the  I'act  that  the  two  sets  of  values  are  in  fairly 
good  agreement.  Also,  the  actual  numbers  are  not  of  that  great  a concern 
since  the  main  purpose  of  this  work  is  to  determine  qualitative  behavior. 

To  completely  detoinnine  the  chemical  kinetics,  the  range  of  tempera- 
ture.s  for  which  the  hydrogen-oxygen  reactions  discus.sed  earlier  are  valid 
must  bo  stated.  .Schott  .and  Kinsey  studied  hydrogen-oxygon  reactions  in  the 
temperature  range  between  1100"  and  26U0®K.  Slightly  broader  limits  were 
chosen  for  this  work,  specifically  IIKD®  to  2700®K.  Again  this  is  a fairly 
arbitrary  choice  but  the  limits  used  .orrespond  roughly  to  values  of  three 
to  nine  times  tho  usual  standard  atmosphere  aml>ient  air  tomperaturo  at  sea 
level.  In  this  tomperaturo  range  thee,  the  .simple  Arrhenius  kinetics  are 
assumed  to  hold  for  the  chemical  process  under  consideration,  below  IUOO“K, 
the  temperature  is  assumed  to  bo  lew  uiough  that  no  appreciable  reaction 
will  occur.  Alwve  27U0®K,  there  is  ai  sumed  to  he  so  much  endotliermic 
dissociation  that  no  appreciable  fraction  of  the  total  delay  time  will  be 
used  up.  This  situation  is  sbown  in  I Ig.  u. 

Next  it.  is  necessary  to  derive  the  specil'ic  form  of  eqa.  8 which  is 
ur.able  v^iith  dl.OUl)  program  ouljtiu  fur  nitlation  and  iU*lay  time  studies. 

Since  dl.UUU  uses  non-dimennion il  val».i*s  of  time  the  .taantity  dt  is  not 
directly  obiainabie.  liuvkever,  llte  qu.intity  dt/t  s *btainable  where "t. 


.'it. 
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has  been  defined  previously.  Thus  onij  can  write 

Xft,  to  ^ 


CIO) ; 


Now  all  that  remains  is  to  substitute  the  expression  forT'^  into  eqn.  10. 
However  is  given  in  terms  of  the  oxygen  concentration  and  the  tempera- 
ture of  a cell  of  the  combustible  mixture.  Thus  expressions  must  first 
be  found  for  these  two  quantities. 

To  determine  the  cell  temperature,  one  may  proceed  as  follows.  Since 
for  a perfect  gas  the  temperature  is  a function  of  the  internal  energy 
only,  we  see  that  the  cell  temperature  can  be  obtained  from  the  cell  energy, 
which  is  a part  of  the  CLOUD  program  output.  As  noted  earlier  CLOUD  yields 
a value  of  L = ® actual  internal  energy  per  unit  mass 

and  and  are  non-dimonsionalizing  reference  quantities.  It  can  also 


be  shown  that  for  a perfect  ga.s 


e s 


P V 


Yo- 1 

We  may  then  substitute  eqn.  11  into  tlie  CLOUD  expression  for  cell  enorgy 
to  obtain, 


Ul) 


E - 


Po  Vo 


to“i 


(12) 


wliich  by  the  perfect  gas  law  miy  be  re-wrltton  as 


c - 

t - 

f<o  ‘o 


va 


whore  is  the  iuiliul  ambient  tempi*'*ature  and  is  the  gas  ooastuut  of 
the  ambient  gas.  Uut  since  we  are  oonsldoriag  vlst  the  ambient  gas  is  the 


10 


'f~-  s 


|)!  ; 


combustible  mixture  we  can  say  that  R = and  eqn.  13  reduces  to 


E = 


T,  Ko- 1 

or  rearranging,  we  obtain  an  expression  for  the  cell  temperature  as 


T=  ET„(!!o-1) 


(14) 


Finally  the  oxygen  concentration,.  [0.,]  must  be  determined.  Since 


S6C*ni0l6S 

A,  the  pre-exponential  factor  \s  given  with  the  units  — , we 


require  to  have  [O2]  in  the  units  so  that  in  eqn.  7 we  will  be  left 

with  units  of  time  only.  The  expression  for  [0^]  is  then  given  by 


" .( Plifir)/ (MW 


(IS) 


whore  ^ and  MW  refer  to  the  density  aud  molecular  weight  of  the  coll  in 
qviestion  and  M.F.  refers  to  the  mole  ''raction  of  O^. 

The  density  of  a coll  is  not  obtainable  directly  from  CLOUD.  However 
the  specific  volume  ratio  v/v^^  is.  Then  because  ^ = 1/v  we  have 


-|-=  s x/(^) 


whoi'c  again  the  subscript  x'efcfs  to  initial  amnient  coiulitions.  Ihen  p is 


given  by 


P = 


Ub) 


Uqn.  16  is  then  substituted  bad  into  eqn.  l!i  to  obtain  the  e.xpre.'ision 


fox'  (O^J , 


(17) 
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Eqns.  14  and  17  may  then  be  subst  Ltuted  back  into  eqn..  7 to  obtain 
an  expression  for  T.» 

n =r7 wT^ r-M — wp(Efl/!?ETo(fc-l)) 

exp(E<./l?£To(Vil)  CIS) 

We  can  now  substitute  intc  eqn.  18  numerical  values  for  the  following 

quantities,  A»6.?4‘l^x40'“ (f  rem  ^IreKlow) 

”w^"“  ('f'r&m  Si'r^Klow) 

MWmjwv«8  ® MWai^  a 

s 4,.^A4*7  y,ter 

Rv00X766  

1.'4* 

Using  those  values  then  eqn.  li>  reduces  to 

Ti^TTOST  » 10-’  (-^)  exp(6?.589eA)  US) 


where has  units  of  seconds,  llucall  also  that  U is  the  non-dimensional 
cell  energy  obtained  from  CLOUU. 

liqn.  19  is  novs  subst itutt^il  back  into  eqn.  10  to  obtain, 

loM»to-’l^xl(«iTw7£T  ' to^ 


~~m tT  '1^ 


where  iite  factor 


1 . 2976  .X  10*^  e”  3698/B) 

L*  I a ■ O' 


i«-\) 


qua  I to  !/T^. 


bqn.  20  is  now  integrutetl  nuuuu'it ally  from  the  nun-dimensional  time 
of  shock  passage  to  the  tiwei%  . la*r  this  calculat  ion,  inpu  20  is  re- 
written  as: 


fU<7T6xio«)e''‘*'”’*'‘'  At  . 1 i,r 


or  more  compactly, 

i/tc 


W 

TtT  r">  -i 

T «iout>  X , r + ” -f,.  ’ 


where  l/l*^  has  been  defined  above.  Note  that  the  quantity  V^^^loUD  units 
of  soc"^  and  that  V'  is  dimensionless. 

The  actual  integration  is  then  performed  using  eqn.  22  or  equivalently 
eqn.  21  on  an  IBM  360  computer.  Tlio  limits  of  integration  arei./i  und 

S 0 

as  shown  above.  However,  it  mu.si  be  noted  that  eqn.  22  was  derived 
with  the  assumption  that  the  reactive  kinetics  are  valid  only  between  the 
temperatures  of  1000®  and  270U“K.  If  the  cell  temperature  drops  below 
1000®K,  the  integration  is  tomj>orarily  halted.  If  the  coll  temperature 
riso.s  above  2700“K,  the  integration  i:  likewise  temporarily  halted.  If  and 
when  the  coll  temperature  retu.'ns  to  the  acceptable  ra.igo,  the  integration 
i.s  then  rosuiacd.  Thus  the  actual  limits  of  integration  aiay  bo  luu'rower 
than  those  given  al'ovo. 

In  I'ig.  7 are  shown  qualitatively  four  diflerent  types  of  temperature 
variation  which  may  bo  observed  for  a given  coil.  In  I'ig.  7a,  tlie  shock 
pusses  the  cell  in  question  and  raises  its  temperature.  However , the 
temperature  never  reaches  the  I000®k  point.  This  cell  will  then  never 
begin  reaction,  itsV'^ulue  wiU  rema  n cero  an^  it  will  have  an  Infinite 
delay  time,  i.e. , it  will  net  ignite,  big.  7b  shoves  a cell  for  whicii  the 
temperature  does  rise  above  IbUU'K  alter  shock  passage.  Once  past  the 
point  tif  highest  tt>mperature,  there  iv  a monotonic  decrease  until  the 
temperature  is  onee  again  below  ihe  ; )U0*^k  cuto'  f limit.  While  tins  cell 
has  a temperature  in  tlie  aceenttiMe  ringe,  iniej'i'aium  of  eqn.  22  along  its 


I ^ 


trajectory  will  proceed.  Thus  this  cell  will  undergo  reaction  and  have 
finite  values  of^Y*  and  delay  tine.  Tills  cell  may  or  may  not  reach  the 
point  of  ignition.  This  will  depend  on  the  initial  dimensional  condi- 
tions and  will  be  discussed  later.  Fig.  7c  shows  a cell  for  which  the 
temperature  rises  above  the  upper  cutoff  limit.  This  cell  will  undergo 
the  same  procedure  as  that  for  the  cell  of  Fig.  except  that  during  the 
time  when  the  temperature  is  above  the  27U0“K  point,  the  integration  of 
eqn.  22  will  be  temporarily  halted,  ihis  integration  will  resume  when  the 
temperature  has  dropped  below  270i)°K  .ind  will  proceed  as  long  as  the  temper- 
ature is  greater  than  1000“K.  All  else  is  the  same.  Finally,  Fig.  7d 
shows  a terapei’aturo-time  plot  for  which  non-ideal  effects  are  pronounced. 

Due  to  the  non-ideal  natui’e  of  the  flow,  second,  third,  etc.  shocks  may 
propagate  past  a cell  after  the  load  shock  has  passed.  The  effect  of  these 
later  shocks  is  to  cause  a tomporatun  pei'turbation  on  the  generally  de- 
creasing temperature  profile.  These  perturbations  raise  the  temperature 
of  the  ceil  but  not  by  as  much  as  the  lead  sliock  raised  the  temperature. 

As  shown  in  the  plot,  a perturbation  lun  occur  in  .such  a way  us  to  t'aise 
the  temperature  above  the  upper  limit  and  so  cause  u temporary  halt  in  the 
integration.  It  can  also  raise  the  temperature  above  the  lower  limit  and 
hence  allow  a longer  integration  for  a cell  fur  which  the  integration 
would  otherwise  hu\c  been  completed. 

With  this  information  in  mind  tl>i>n,  eqn.  2.  is  integrated  along  the 
neco.ssury  particle  or  cell  tv.jei‘tori*'S  to  obtain  the  data  for  this  investi- 
gation. 


IV.  fl'li'i  CASli 


l‘o  determine  the  validity  of  vlu*  uumorical  integrating  seheme,  a test 
case  was  nm  witli  the  Huang  and  HIuhi  .'(Hi  atm  bursun,,  sphere  run  described 
earlier.  The  numerical  iiuegr.at ing  .‘heme  was  ‘un  on  a nuisiber  of  cells 
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for  this  sphere  burst.  iVo  of  these  i.Jll?  are  discussed  here,  namely 
cell  130  and  cell  135.  For  these  two  cells,  the  quantity  1/T^  was  also 
obtained  at  various  times  during  the  integration  period.  Plots  were 


then  made  of  vs  t/t^  for  these  colls.  These  are  shown  in  Figs.  8 


and  9.  As  eqn.  22  can  also  be  integrated  graphically  by  determining  the 


area  under  the  1/t^  vs  t/t^  graph  this  was  also  done  and  used  as  a check 


on  the  numerical  scheme.  In  Figs.  8 and  9 the  cross-ltatched  areas  represent 
half  time  increments  which  were  included  to  make  the  graphical  integration 
correspond  exactly  to  the  numerical  integration.  The  area  under  tliese 
ctu'ves  was  tlien  deteimiined  by  a plani'i>eter. 


For  cell  15U,  the  numerical  scheii'C  yielded  a value  {^lOUD 


4-1  4 

3.379  X 10  .sec  . The  planimtter  integration  yielded  a value  of  3,376  x 10 


-1 


sec  . The  percent  error  can  tlien  be  determined  from. 


% error  = 


■X 100 


(23) 


For  cell  150,  eqn.  23  yields  a value  I'f  .09‘»  error. 

F'or  coll  13.5,  tito  numerical  sci'cme  yielded  a value  of 

4 “ I 4-1 

1.622  X 10  sec  whox'eas  the  planimeter  yielded  a value  of  1.62  x 10  sec  . 

lU|n.  23  then  gives  u value  of  on  or  for  cell  135. 

Tho.se  values  of  pereontago  error  are  certainly  within  an  acceptable 

error  marglti  and  so  thi.s  is  ta^o'  as  « vidence  that  tlie  numerical  integration 

is  returning  valid  output.  / 

Also  shown  for  ilUisU’ative  purjose.s  is  Fig.  10,  which  is  a three- 

dimensional  plot  of  Tfll  uUP  ^ .•ell  muiher  vs  obtained 

with  the  numerical  integration.  Note  that  for  all  the  cells  plotted,  a 

large  fraction  of  the  final  valiu  of  obtained  very  quickly,  even 

for  the  first  four  cells  in  which  tlu<  integration  proceihis  for  a long  time. 

Note  also  maN  nuain  a value  greatv  than  unity  even  though  ^ 


.yy ^ . 


H : 

i 
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itself  is  limited  by  this  value.  This  fact  becomes  important  for  ex- 
plosion loci  determination  and  will  be  discussed  below. 

It  is  desired  to  determine  three  things  for  the  liuang  and  Chou  200 
atm  case;  the  delay  time  to  ignition  for  various  cells,  explosion  loci 
and  minimum  energy  requix aments  for  ignition. 

First  of  all,  the  value  of  was  obtained  for  a number  of  cells 

for  this  I'un.  These  values  are  plotted  vs  cell  number  on  Fig.  11.  Notice 
again  greater  chan  unity,  l-qn.  22  can  be  re-written 

simply  as; 

Now  for  a coll  to  ignite,  its  value  ol'Y' west  bo  equal  to  unity.  Then 
oqn.  22  reduces  to: 


Taoui>*  t« 


or 


U4) 


aeui> 


Thus  knowing  the  value  V*ci OUU'  value  of  t^,  which  ha.s  been  used  to 
non-dimonsionaliao  tho  time  va-iahle  can  be  dettnuined.  Then,  knowing 
the  value  of  t/t^  at  which  Vi''0UI>  causes  V' 

for  the  cell  to  equal  unity,  one  can  easily  deunnine  the  real  time 
^'Dl’LAY  ignite  iTom: 

toei«  --  to  (25) 

Suppose  now  that  one  looks  at  th>  ma.tiinum  value  li* 

This  is  equal  to  .\  lU^’  see'^  and  occurs  tor  cell  82.  Hubstituting 

this  value  into  eqn.  24  one  finds  th.u  « S.u'.i'.'  x it)''  sec.  From  the 
computer  output,  It  is  known  that  cell  82  sums  to  a \alue  '■’i'lfpiQyp 
= l.‘Jt)2  X U)^’  sec"^  at  non-dimensional  t/t^^  I Thus  from  eqn.  21>, 


■ill 


for  coll  82 


"toUAV  - (1.15  S'!  (.5.0 77  >10*'’  S€C.) 
^S.ni  >10-''  iec. 


If  the  value  o^VcLOUD  value  of  t^  just  determined 

are  entered  back  into  eqn.  22,  they  will  of  course  yield  a value  ofT=i 
for  cell  82,  However,  if  the  value  of  ether  cell  in  the 

Huang  and  Chou  run  i.s  entered  into  eqn.  22  with  the  given  t^,  it  will 
yield  a value  V<  1.  Thus  onl)  coll  »2  will  undergo  ignition. 

Suppose  however,  that  one  considers  a value  which  is  less 

than  the  maximum  value  'VGlOUD'  *1K^L01)D  ^ ^ > ''^^ich 

is  the  value  obtained  for  cell  85,  Then  from  eqii.  24,  we  obtain  t^ 

-7  -7 

= 5,591  X 10  sec.  From  eqn.  25  we  ^.ot  tor  cell  85  » 6,119  x 10  sec, 

where  the  value  of  t/t^  for  "'1'''^^' ^(;i0Ul)  completely  summed  was 
obtained  from  the  computer  output.  Again,  entering  the  value 
for  cell  83  with  the  just  dutormimd  hack  into  eqn.  22  gives  Y.  1,  as 
o.\pect'.ed.  Suppose  that  one  now  looks  at  ceils  80,  81,  and  82.  Those 
cells  have  values  wh.ch  are  greater  than  the  value  obtained  for 

coll  83.  So  when  those  values  are  «n)erod  back  into  eqn,  22  with  the 
given  above,  they  will  yield  values  vhicl»  are  greater  than  unity.  This 
muun.s  that  in  addition  to  cell  83,  tiu'se  other  throe  cells  will  also 
ignite.  Of  course,'^/'can  not  physica.ly  exceed  unity.  What  is  done  to 

correct  this  is  to  look  at  the  value  »*f  i/t  for  whlc.'i  each  of  the.se  cells 

* 

attained  the  t’^ltje,  jup  * 1.851*  Id 'sec'  . This  is  of  course  the 
value  ViTOUH  summation  beyond  tliis  value  is 

then  Ignored.  Now  all  of  tl^es^>  cells  will  ytelc'V'*  1 In  eqn.  22. 

One  can  of  cowr.se  consider  vaUn,  ■>  ot  "^Cl.ui'b  wh i ch  are  still  1 ower , 
e.g.  ^ l.elib  ,\  Irt^scc  the  ’uhic  oht.iined  lor  cell  85.  Now,  all 

the  cells  from  79  to  85  v»ill  i.;n*U-.  llicir  various  ignition  puiametors 


may  be  determined  as  above. 

Since  each  value  1^'L(JU£  uniquely  determines  a value  of  we 

can  now  work  back  to  obtain  dimensional  initial  conditions  for  the  bursting 

14 

sphere.  This  is  done  as  follows.  Adamczyk  shows  that  there  is  a 
relation  between  t^  and  r^^,  the  initial  sphere  size.  This  is  given  by: 


fo  '/T 


where  a^  refers  to  the  ambient  undisturbed  sound  speed  and  Yj  is  as  defined 
earlier.  For  a perfect  gas,  the  quantity  a^  can  be  replaced  by, 

RoTo  (27) 

Thus 

, r./^ 

However,  since  this  investigation  assumes  air  on  both  sides  of  the  contact 
surface  (for  flow  field  determination)  wo  have  that^^  3 e 1.4  and 

■t  - 

• yiTir’ 

which  can  bo  ro-writton  as 

r.»to/R7fr  (28) 

where  is  determined  from  ignition  uualysi.s  as  above,  ° 


» 6.SSS  X 10 


and  « 208**K.  liqn.  28  yield.s  the  initial  radius, 


r^  of  the  bursting  sphere. 


Uno  can  then  use  this  value  of  r in  eqn,  4, 


to  obtain  the  quantity  H^,  which  as  discussed  earlier,  is  the  amount  of 
energy  which  must  be  added  to  yield  the  initial  sphere  burst  conditions. 
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Table  2 displays  the  ignition  data  discussed  above.  For  three  values 

°^X:L0UD‘  values  of  and  shown  for  cells  79 

through  85.  Also  included  arc  tae  values  of  t , r and  E,  for  each  case. 

* o o 1 

TABLE  2 


c 

’^cloud“  ^ 

.962x10^  sec*^ 

'^CLOUD'  ^ 

.855x10^  sec’^ 

'^CLOUD' 

656x10^  sec'^ 

E 

t 5.097x10*^  sec 

0 

= 5.391x10"^  sec 

t = 6*059x10’"^  sec 

L 

r = 1.491x10'^  cm 
0 

r = i, 577x10"^  cm 
0 

= 1.766x10'^  cm 

L 

Ej  a 6.992x10'^  joules 

Ej  = 8.273x11)'^  joul  .'S 

Ej  = 1.163x10'*^  joules 

t /t  1 

*DUAY  t*®'’ 

79 

00 

00 

to 

«) 

.290 

1.75x10’^ 

80 

QO 

« 

.*.75 

•1.178xlU‘‘ 

.230 

1.389x10’^ 

81 

0» 

00 

.540 

1,835x10' ' 

.201 

1.214x10’^ 

82 

I.ISS 

5.887x10*^ 

.500 

1.617x10'’ 

.195 

1.178x10'^ 

85 

a* 

0l> 

1.I35 

('.119x10''^ 

,227 

1.571x10'^ 

84 

m 

.295 

1.782x10'^ 

85 

eo 

30 

ili 

1.115 

6.733x10'^ 

I 


t'ig.  12  shows  a plot  of  thu  dimomionui  igiution  doiay  tiutos  vs  coll 
number  for  the  given  values  of  In  the  table  above.  N‘oto  that  cell  82 
has  the  shortest  delay  time  for  all  thiee  energy  level.s.  This  makes 
intuitive  sense  because  coll  82  had  the  highest  value  of However, 
it  may  occur  that  another  cell  will  siu,  to  a given  value  before 

coll  82  does,  even  though  cell  :i2  eventually  sttms  to  the  highest  value. 

In  such  a case,  U that  value  chosen  as  a reference  for  an 

ignition  analysis,  veil  82  wilt  nut  hu'o  the  dh'Uiest  delay  time.  Thus  one 
eaii  not  say  in  general  that  ceil  82  wt  l alway.s  first. 


’tl^l  A|AV|-.|  y-- 
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It  is  also  seen  from  Fig.  12  that  as  the  energy  added  to  the  bursting 
sphere,  increases,  the  delay  time  for  all  the  cells  decreases.  Again, 
this  makes  intuitive  sense  because  one  would  expect  a quicker  reaction  for 
a larger  amount  of  energy  in  the  source  area. 

Fig.  13  shows  a pictorial  view  of  the  Huang  and  Chou  200  atm  flow  field 
generated  by  the  CLOUD  program.  On  this  graph,  the  axes  are  non-dimensional 
distance  r/r^  vs  non-dimensional  time  t/t^.  The  particle  paths  or  cell  tra- 
jectories are  shown.  The  lead  shock  is  very  evident  in  the  lower  left  cor- 
ner. Initially,  the  cells  travel  vertically,  i,e.,  they  cover  no  distance 
with  time  and  are  stationary.  VIhen  the  load  shock  arrives,  it  bends  the 
trajectories  towards  the  horizontal,  i.o.,  it  accelerates  the  particles 
and  hence  causes  a change  in  distance  with  time.  As  cun  be  seen,  the  part- 
icles eventually  decelerate  and  tend  bi*ck  toward  a vortical  inclination. 

As  mentioned  earlier,  duo  to  the  non-ideal  nature  of  the  flow,  secondary 
wave  systems  will  occur.  These  can  be  seen  as  ripples  superimposed  on  the 
general  behavior  of  the  coll  trajectories. 

The  blockod^off  portion  of  Fig.  13  is  the  area  of  the  flow  field  which 
is  of  interest  for  ignition  calculations.  This  region  is  magnified  and 
displayed  in  Fig.  14.  Recall  that  in  Table  I the  quantity  was 

tabulated  against  cell  number  for  various  energy  levels,  For  each  of 
these  energy  levels,  one  can  plot  the  |*oints  at  which  Ignition  occurs  for 
cells  in  question  by  sim|)ty  loctiting  the  proper  trajectory  and  travel littg 
along  it  until  the  required  value  of  on  the  vertical  axis  is  reached. 
These  iMDints  can  then  be  connec:ed  to  form  what  will  he  called  the  explo- 
sion loci.  These  loci  are  plotted  for  two  energy  levels,  namely 
a g,273  X 10~^  joules  and  F|  « l.lb.i  x lU"*'  Unties  in  Fig.  14. 

N'ote  that  the  explosion  locus  for  the  lower  value  of  the  energy  Is 
higher  up  the  cell  trajectories.  This  corres}Hmds  to  the  fact  that  a lower 


V.  BURSTING  SPHiiRti  RESULTS 

The  above  procedure  described  for  the  Uuun;',  :iud  (.:hi>u  case  has  also  boon 
applied  to  a number  of  other  bui’sting  ; phere  runs.  Those  sphere  Iniv.sts  are 
generated  by  the  CLOUD  pi’ogram  iind  are  characteiM  :ed  by  three  parumetevs: 
the  pressure  ratio  P^/Py.  tl»o  tomporature  rat  io  I aiul  the  values  of  ^ 
for  the  source  region  and  aiiihiont  region,  called  ^ , and  respectively. 

For  all  the  flow  field.s  goneratod  by  the  CLOUD  pj'ogriim.  the  source  and 
ambient  regions  are  assumed  to  contain  air.  Thus  Yj  =1.4. 

It  is  first  necessary  to  determine  which  combinations  of  these  purameter.s 
will  load  to  flow  fields  in  which  init'Ution  may  occur.  Tliis  criterion  cun 
bo  determined  analytically  in  the  following  inaniur. 

Wo  have  sot  the  minimum  temperature  for  reaction  at  1000*K.  Therefore 


. energy  will  result  in  a longer  time  for  reaction,  as  was  discussed  earlier. 

-4 

The  explosion  locus  for  the  energy  level  = 6.992  x 10  joules  is  not 
shown  as  it  is  only  a single  point. 

By  this  tiur,  it  should  be  evident  that  there  is  something  special  about 
-4 

the  = 6.992  x 10  joules  energy  level,  because  it  allows  only  one 
particle  in  the  entire  flow  field  to  ignite.  What  this  means  is  that  this 
is  the  minimum  amount  of  energy  which  will  result  in  an  ignition  in  the 
Huang  and  Chou  200  atm  flow  field.  If  this  energy  were  decreased  slightly, 
it  would  cause  every  particle  in  the  flow  field  to  have  an  infinite  delay 
time  which  means  of  course  that  there  would  be  no  ignition  of  any  kind. 

Thus  this  value  of  represents  a cutoff  between  ignition  and  no  ignition 
for  this  particular  flow  field.  The  values  of  the  lowest  energies  for 
ignition  will  be  obtained  for  other  flow  fields  in  the  following  sections 
and  the, general  behavior  of  thei-o  quantities  will  be  di.scussed. 
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the  flow  field  must  show  temperatures  at  least  this  high  behind  the  incident 
shock.  Now  for  a bursting  sphere,  the  shock  will  have  maximum  strength  at 
the  instant  of  burst.  From  this  point  on  it  will  decay,  eventually  to 
approach  an  acoustic-type  behav;.or.  Thus  at  the  moment  of  burst,  the  shock 
must  be  capable  of  raising  the  temperature  of  the  mixture  to  1000°K.  One 
can  now  use  the  Rankine-Hugoniot  relation  for  temperature  from  Liepmann 
and  Roshko.  , Ps 

T,  _ Ps  , 

TTnor 

**  TTj-  p, 

where  is  the  temperature  ratio  across  the  shock,  Pg/P^  is  the 

pressure  ratio  across  the  shock  and  ^ replaces  and  since  it  was 
assumed  that  ? ^ = X = 1.4.  If  = 298®K  and  at  the  moment  of  burst 

wo  must  have  1000°K  behind  the  shock  tlion, 

Ts  . aOOO*j< 

■X"  -STS*'' 

Wo  can  then  substitute  this  value  of  T»/T^  and  U = 1.4  back  into  oqn.  29. 
After  .some  rearrangement  it  is  found  that 


|■^l+.134a{-|•)-3.355T  »0 


This  i.s  u quadratic  in  the  quantity  H way  bo  solved  to  yield  the 


values 


•^=14,3«7%  ,-0.^i'S36 


whoro  the  negative  solution  is  not  physically  meaningful.  We  now  have  the 
pressure  ratio  across  the  shock  at  the  instant  of  hurst  which  is  necessary 
to  produce  a lOOU'^K  toiuporature  Recall  eqn.  o 
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'a»)(P5/fe-i) 


+ (5  + 


where  y has  now  replaced  Yj^  and  y^.  The  quantity  C^q/^2^  expressed 


iQ. 

V 1ft  R 

Again  because  we  are  considering  air  for  the  flow  field  calculations,  the 
specific  heat  ratio  and  gas  constant  ratio  drop  out  and  we  have 

-3^  c (Jsjf' 
fiA  \ J 


do  ^ / T>  yv^ 
aa  ' \ % I 


Substituting  eqn.  30  back  into  eqn.  6 we  obtain 

r /v  <tVf!  /r  \"*/a  / 


^ - aF  1 _ 

15 " Po  ^ 


which,  after  substituting  in  tlie  values  determined  for  y and  reduces 


w 


T»v'^r^ 


liqn.  31  now  relates  the  initial  sphere  pressure  ratio  to  the  initial 
sphere  toroporaturo  ratio,  while  satisfying  the  constraint  of  a 1000“K 
temperature  behind  the’ lead  shock  at  the  instant  of  hurst. 

Suppose  tlie  minimum  tompu;.*ature  liufined  for  the  Arrhenius  reaction 
kinetics  had  boon  1200'*K  instead  of  in00'’K.  Although  this  criterion  is  nut 
made  use  of  in  this  investigation,  it  can  he  mentioned  for  comparison 
purposes.  If  this  were  the  ca«e  however,  then  an  equation  similar  to 
eqn.  31  could  be  derived  using  the  sane  logic.  This  equation  would  he, 


[l-.6W7 


Eqns.  31  and  32  are  plotted  on  Fig.  15  for  comparison  purposes.  For  a 
1000®K  minimum  temperature  it  can  be  seen  that  only  points  which  are  to 
the  right  and  above  the  1000®K  curve  will  provide  the  proper  initial  sphere 
conditions  for  an  ignition.  No  set  of  initial  sphere  conditions  below  and 
to  the  left  of  the  curve  will  produce  an  ignition  in  the  flow  field.  A 
corresponding  argument  can  be  made  for  the  1200“K  curve. 

Note  that  the  curve  for  a 1200®K  minimum  temperature  lies  higher  on 
the  graph  than  the  curve  for  1000“K.  This  makes  sense  because  one  would 
expect  that  it  would  require  larger  initial  values  of  sphere  temperature 
and  pressure  ratios  to  cause  a higher  temperature  in  the  flow  field.  Note 
also  that  for  either  curve  as  the  temperature  ratio  decreases  to  zero  the 
pressure  ratio  rises  very  rapidly.  Likewise,  the  roquii'ed  temperature  ratio 
increases  markedly  as  the  presiiure  ratio  decreases.  Also  shown  on  this  graph 
is  the  curve  of  temperature  ratio  equc*l  to  pressure  ratio. 

For  this  investigation  a number  of  bursting  sphere  runs  were  made  to 
study  the  ignition  behavior  in  the  flow  field,  Ihe  initial  conditions  for 
all  of  those  runs  lie  above  and  to  tho  right  of  the  lOOO^K  curve.  The 
intersection  of  the  lOOU'^K  curvo  and  tho  cut've  of  temperature  ratio  equal 
to  pressure  ratio  was  chosen  as  a .staj*ting  point.  This  intersection  cor- 
responds to  tho  initial  conditions: 

-1-1 30  -i-sJO 


As  thLs  set  of  initial  conditions  lien  on  the  cutoff  curve,  it  was 

not  expected  to  yield  any  uctu.il  finite  ignition  data  with  the  mauerical 
routine.  Thus  to  provide  the  'nutueric.tl  integration  with  a flow  field 
slightly  above  the  cutoff  flow  field,  the  above  initial  conditions  were 
increased  by  ten  percent.  Thu* 
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= 33 


?o  To 

Then  with  this  set  of  more  reasonable  initial  conditions  as  a 
starting  point,  eight  bursting  sphere  flow  fields  were  generated.  Four 
of  these  flow  fields  successively  doubled  the  initial  pressure  ratio  while 
holding  the  initial  temperature  ratio  constant.  The  other  four  flow  fields 


successively  doubled  both  the  initial  pressure  and  temperature  ratios. 
These  flow  fields,  or  equivalently  computer  runs,  as  they  will  henceforth 
be  called  are  displayed  in  Tables  3a  and  3b. 


i"- 


FLOW  FIELD/RUN PRESSURE  RATIO  THMPERATURF.  RATIO 


THEORETICAL  CUTOFF 

30 

30 

A)  10%  ABOVE  CUTOFF 

33 

35 

1 

66 

53 

2 

132 

53 

3 

264 

53 

4 

528 

35 

(aj 


FLOW  FIELD/RUN 

PRESSURE  l-ATIO 

■TEMPERATURE  IbVTIO 

THEORETICAL  CUTOFF' 

.30 

50 

A)  10%  ABOVE  CUTOFF 

53 

33 

5 

66 

b6 

6 

152 

132 

7 

264 

- 264 

8 

.528 

528 

(b) 


FLOW  FIELO/RUN 

PUUS.SUlUi  RATIO 

ti:mpeiu‘ture  ratio 

THEORETICAL  CUTOFF 

60 

10 

B)  10%  ABOVE  CUTOFF 

66 

n 

0 

152 

n 

10 

264 

11 

11 

S28 

n 

i's) 
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It  was  also  desired  to  study  ignition  behavior  in  another  portion  of 
the  admissible  region  of  Fig.  15.  For  this  purpose  a new  point  on  the  1000®K 
cutoff  curve  was  chosen  as  a starting  point  for  another  series  of  runs. 

The  initial  conditions  at  this  cutoff  )>oint  were 


A.- 
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As  before  these  values  were  increased  by  ten  percent  to  yield  a better 
starting  point  for  numerical  integration.  Thus, 


Jk_. 


= 1J 
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Then  with  this  new  set  of  initial  conditions  as  a starting  point, 
three  more  bursting  sphere  runs  were  made  in  vjhich  the  pressure  ratio 
was  successively  doubled  while  the  temperature  ratio  was  held  constant. 

These  runs  ui'o  also  shown  in  Table  3c. 

As  was  done  fur  the  Huang  and  Chou  200  atm  case  in  Fig.  11,  tho  values 
of  l^Louo  number  are  plotted  lor  runs  1 through  11  in  Figs,  lb,  17, 

and  18.  Two  things  can  be  observed  in  Figs.  16,  17  and  18.  First  for  a 
given  sot  of  runs,  tho  maximum  value  ul  "^L0UI»  utcreuse.s  us  tho  inin  number 
ineroasos.  Then  since  the  initial  pres.suro  ratio  increases  as  the  run 
number  increases,  we  conclude  that  the  luaximum  value  inct'easos 

us  the  initial  pressure  ratio  Increases. 

We  would  now  like  to  relate  this  pressure  dependence  of  the  maximunt 
value  an  energy  dependence.  This  can  be  done  with  the  fol- 

lowing derivation.  Lot  the  quantity  1*  ,/P  » x.  Then  we  may  write 

* V 


P, » P.  X 


and  we  derive: 
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But  X =5  ^ 
0 


Thus  we  can  see  that  an  inci'easo  in  the  initial  pressure  ratio  yields 
a raonotonic  increase  in  the  expression  Now  trom  eqn.  4 we  know 

that  the  initial  sphere  energy,  is  directly  proportional  to  the  quantity 


P^  - P^.  Thus  wo  may  finally  conclude  that  the  maximum  value  of  ij; 

*L  0 ' 
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increases  as  the  initial  sphere  energy  increases.  Tliis  is  a logical 

conclusion  because  it  has  been  shown  that  a largo  value  of  yields 

a small  value  of  t and  henco  small  values  of  dimensional  delay  time.  Thus 
0 

wo  see  that  as  the  energy  increases,  the  value  of  the  maximum 
will  increase  and  corresponding  to  this,  tiic  dimensional  delay  times  will 
docroaso.  This  is  exactly  wiiat  one  would  expect.  It  should  bo  noted  how- 
ever that  there  may  bo  slight  perturbations  on  this  general  trend  due  to 
tlto  non-steady  and  iton-idoul  nature  of  the  flow  fields  under  consideration. 

It  is  interesting  to  note  that  the  quantity 
has  been  studied  previously  by  Strehluw  and  Uicker  and  Adamcsyk^'*.  This 
quantity  is  called  the  dimensionless  energy  density  for  the  bursting 
sphere  and  is  given  the  symliol  q.  Thus 

P 

s — - 1 i34) 

The  quantity  q represents  the  lueunt  of  energy  added  to  the  sphere 
normalised  by  the  amount  of  energy  tint  would  have  heeit  in  an  equivalent 
volume  but  unpressurized  sphere.  Proa  eqn.  t ^ is  given  by, 


Cl  "Hr  Tf,.i  • 

- fP*i"  V 
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The  energy  for  the  equivalent  volume,  unpressuri zed  sphere  will  be  given 
by  the  quantity  PoV^/(Yj  - 1).  Then  we  have 

P> " M -1 


■ P.V. 


0 Vo 


■ ^^TT^'^Okwc.  P, 
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or 

The  quantity  q will  be  used  later  in  thijj  investigation  when  delay  times 
and  minimum  energy  requirements  for  ignition  are  studied,  because  it  has 
a closer  physical  relation  to  the  energy  than  does  the  actual  initial  sphere 
pressure  ratio. 

The  second  noteworthy  fact,  in  ri'gayd  to  l*igs.  lb,  17  and  18,  is  that 
the  cell  which  has  the  largest  value  of  '•'an  not  ho  determined  prior 

to  the  actual  numerical  integi'ution.  Thl^  i»  because  the  cell  for  maximum 
^’10110  **  complicated  function  of  the  flow  field  variables,  which 

are  themselves  complicated  functions  oi' initial  bursting  sphoi'e  conditions. 
From  oqii.  21  one  may  note  that  for  a eell  in  a given  run, 
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(i) 


Thus  is  a complicated  fanction  of  the  time  histories  of  the  dimension- 

less cell  energies  and  specific  volumes, 

One  can  however  state  that  the  mtxii^niii  va.ue  of  does  tiot  in 

general  occur  at  or  around  the  cent aet  cell:«  SO  and  SI).  For  runs 

1,  2 ujid  5 the  coll  which  has  the  lar;os|  valuv  of  iscoll  51  which 

is  adjacent  to  the  contact  surface.  iv‘hh|  is  irohahly  luippeniiig  tn  these 
runs  l.s  that  as  one  progres?»es  out  fivw)  the  contact  surface,  the  coll  energy 
F.  decays  rapidly.  Then  the  tetia 
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gets  very  ssiall  and  results  in  a small  value  of  ♦(.^jup.  This  behavior  ts 

not  evident  for  other,  higher  initial  energy  runs. 

Finally  in  regard  to  Figs.  16.  17  and  18  one  fact  must  be  noted.  For 
a few  of  the  cells  in  runs  7 and  8 the  numerical  integration  could  not 
proceed  to  completion.  This  was  due  to  the  fact  that  in  these  runs  the  cell 
energies  became  so  high  that  the  temperatures  did  not  decay  below  the  1000»K 
cutoff  value  by  the  time  the  CIOUO  program  reached  the  limits  of  its 
array  dimensions.  For  some  analyses,  this  would  not  lead  to  large  errors 
because  as  was  evident  from  Fig.  10,  colls  attain  most  of  their  final 

values  of  *(;louu 

field  is  seen  to  contribute  only  a small  amount  to  However  for 

the  delay  time  and  minimum  ignition  energy  analysis  discussed  later  the 
above  incomplete  numerical  integration  could  lead  to  some  er.-onoous  results. 

This  will  bo  discussoU  detail  later. 

Further  data  reduction  wua  undertaken  tor  runs  I through  4.  Tables 
4,5.6  and  7 shew  data  for  these  runs  as  was  sl>ow«  for  the  Uuoitg  and  Chou  200  atm 
case  in  Table  2.  Those  results  are  then  plotted  in  Figs.  10»  20,  21  and  22 

as  was  also  done  for  the  Huang  and  Chou  case. 

Figs.  10,  20,  21  and  22  show  an  interostini;  range  of  delay  time 
behavior.  For  runs  I and  d.  Figs.  IF  aiul  dd  show  only  half-curves.  This  is 
of  course  duo  to  the  foot  that  in  ilu  ie  runs  the  cells  which  had  the  largest 
values  of  wore  adjaceni  to  ihv  contact  surface.  In  Fig.  19  U cun  also 

he  seen  tliat  cell  ill.  adiatrtit  to  the  contact  surface,  ignites  at  virtually 
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t = 5.793x10  ^ sec  t - 7.776x10  ^ sec  t = 8.482x10  ^ sec 
0 0 0 


r = 1.694x10'^  cm  r = 2.27  IxiO'^  cm 
0 0 
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0 
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'^CLOUD"  2.846x10^  sec'^ 
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0 
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0 
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0 
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• 
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.345 
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.575 
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.500 
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.355 

1.177x10’" 

.275 

1.009x10’" 

.350 

1.230x10*' 

.290 

1.064x10’^ 

.;75 
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t = 2.290x10'^  sec  t = 2.641x10'" 

0 0 
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85  1.175  2.691x10' 

84  ^ A> 
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86  * *** 
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.‘iOS  1.070x10”*' 
.475  1.254x10”^ 
.625  1.651x10"^ 
..»55  9,580x10  ^ 
.595  1.045x10“' 
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.555  1.415x10“^ 


1.  15 


». 945x1  o' 


'^CLOUD"  5*724^xTo^"^'^ 

= 2.685x10'^  sec 

r = 7,853x10"*^  cm 
» 

Ej  = 2.708x10'*^  joules 

— 

^DELAY 

2.565 

6.550x10'^ 

.585 

1.054x10'^ 

.450 
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.555 
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1.155x10"^ 
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1.302x10"^ 

. 605 
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the  same  time  for  two  different  initial  sphere  energies.  Figs.  20,  21 
and  22  show  some  ripples  in  tht  curves  which  are  most  likely  due  to  the 
passage  of  the  secondary  wave  systems  and  general  non- ideal  nature  of  the 
flow  fields.  Fig.  22  especially  shows  some  ratlier  striking  behavior 
around  the  vicinity  of  cell  80.  where  there  is  a secondary  extremum  in  the 
curve  for  each  of  the  two  larger  energy  levels.  This  type  of  behavior  would 
be  almost,  if  not  impossible  tc  deduce  analytically,  whereas  the  numerical 


procedure  shows  the  results  quite  well.  * ^ 

To  compare  these  runs  with  each  other  in  terras  of  delay  time  to 
ignition,  Fig,  23  has  been  obtained.  On  this  graph  are  plotted  the 
delay  times  which  coi'respond  t(>  the  minimum  energy  for  ignition  level  for 
each  given  run.  These  data  are  plotted  vs,  the  non-dimensional  energy 
density,  q v«rhich  is  used  to  churactorico  each  run.  Huns  1,  2,  3 and  4 
which  have  just  been  investigated  yield  points  on  the  eur>^e  labelled 
T,/T  a 53.  Also  shown  on  tnis  curve  is  a data  point  for  run  A.  Included 

on  this  plot  also  are  two  other  curves . The  bottom  curve  corresponds  to 

mis  4,  3 and  t>  for  which  the  lemperature  ratio  initially  in  the  bursting 
sphere  equalled  the  initial  pressure  ratio.  Huns  7 and  S have  not  been 
reduced  for  delay  time  data  due  to  tin*  incomplete  numerical  integration  of 

%L0Ui>  riuis.  the  top  curve  corresponds  to  runs  U,  ‘d,  lb  and  U. 

In  these  runs  the  initial  bursting  sphere  temperature  ratio  was  held 
constant  at  a value  of  T,/l  « 11. 

Looking  at  each  curve  on  2i  alone,  it  is  obvious  thvU  tUv\ 
time  is  a function  of  the  initiai  Inn  .ting  sphere  pressure 
since  a change  in  the  bursting  sphert  leapeiMiure  ratio  causes  a switch  to 
anotheir  curve  aliugether,  it  is  aiso  seen  that  the  delay  time  is  a function 
of  this  initial  temjieratuiv  ratio,  i ;*oa  observing  the  relative  positions 
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of  the  various  curves  one  could  conclude  that  the  relation  between  initial 
sphere  temperature  ratio  and  corresponding  delay  time  is  highly  non-linear. 

The  physical  reason  for  the  dependence  of  the  delay  time  on  the 
temperature  ratio  can  be  shown  as  follows.  We  know  that  the  energy  added 
to  the  sphere  is  a function  only  of  the  pressure  and  not  of  the  sphere 
temperature.  Thus  the  temperature  does  not  contribute  to  the  bursting 
sphere  energy.  Nevertheless  the  temperature  does  influence  the  flow 
field.  At  the  instant  of  burst  a travelling  wave  system  is  set  up.  The 
lead  shock  propagates  outward  from  the  contact  surface  while  a rarefaction 
fan  propagates  inward  toward  the  center  of  the  sphere.  This  rarefaction 
fan,  being  an  isentropic,  acoustic  wave  must  propagate  at  the  local  sound 
speed,  which  with  the  assumption  of  a perfect  gas  is  a direct  function  of 
the  temperature  in  the  bursting  sphere.  The  energy  available  from  the  high 
pressure  sphere  cannot  be  utilized  to  drive  the  lead  shock  until  the 
rarefaction  fan  has  passed  through  and  hence  released  it  to  the  flow  field. 
Thus  the  rapidity  with  which  this  energy  is  made  available  to  the  flow  field 
is  dependent  on  the  sound  speed  in  the  higli  pressure  gas,  which  is  then  in 
turn  related  to  the  temperatui'o  or  temperature  ratio  in  the  initial  Imr.sting 
sphere.  Thus  changes  in  the  s))hore  temperature  ratio,  while  not  changing 
the  bursting  sphere  energy,  do  alter  ihe  flow  field  enough  to  have  a definite 
effect  on  the  delay  times  for  various  runs. 

To  provide  a qualitative  picture  of  the  ignition  behavior  in  the  flow 
field  of  a bursting  sphere  explosion  loci  have  boon  plotted  for  runs  1,  2 
and  3.  These  are  shown  in  I’igs.  24,  .IS,  2u,  27,  28,  and  21).  Vov  each  run, 
the  entire  flow  field  is  depicted  fiist,  with  tl>e  are.»  of  interest  for 
ignition  studios  blocked  off.  This  blocked  off  area  is  then  magnified  to 
provide  a clearer  view  of  the  detailed  beliavior  being  studied.  U.\ploslon 


loci  for  the  two  higher  values  of  for  each  run  are  plotted  in  the  magni- 
fied graphs. 

In  Figs.  24,  26  and  28  note  the  increasing  strength  of  the  lead  shock 
as  the  run  number  and  hence  initial  sphere  pressure  ratio  increases.  This 
vividly  illustrates  the  changing  flow  field  which  results  when  initial 
sphere  conditions  are  varied.  Secondary  shocks  are  also  visible  on  Figs.  26 
and  28  for  cases  2 and  3. 

In  the  magnified  plots,  Figs.  25,  27  and  29  it  should  first  be  noted 
that  the  rather  irregular  termination  of  the  particle  paths  in  the  upper 
right  parts  of  the  graphs  is  not  physically  realistic  but  results  from 
plotter  limitations.  Figs.  25  and  27  again  show  only  one  half  of  the 
explosion  loci.  As  before  thi.s  is  due  to  the  fact  that  the  runs  for  these 
plots  have  their  maximum  value  of  for  the  cell  adjacent  to  the 

contact  surface.  Fig.  29  for  run  3 shows  the  full  expected  explosion  loci. 

It  is  interesting  to  note  that  the  explosion  loci  are  nearly  parallel 
to  the  actual  particle  trajectories  for  these  runs.  In  other  words,  each 
cell  has  a distinct  non>dimon.s  .onal  time  at  which  it  will  ignite.  One  might 
have  expected  the  explosion  loci  to  lie  horizontally.  If  this  were  the  case 
then  it  would  imply  that  a number  of  cells  would  ignite  simultaneously. 
However,  this  typo  of  behavior  has  not  been  observed  for  flow  fields  studied 
in  this  investigation. 

Finally  we  would  like  to  itudy  the  minimujn  energy  for  ignition  behavior 
for  the  bursting  sphere  cases.  Hecal ' for  the  liuang  and  Chou  200  atm  case 
that  there  was  one  value  of  for  which  only  one  cell  in  the  flow  field, 
thni  one  with  the  largest  valuj  of  would  ignite.  This  was  culled 

the  minimum  energy  fur  ignitiot  as  energy  less  than  this  amount  would  cause 
all  cells  to  have  an  infinite  ielay  t tme.  A corresponding  energy  level  can 
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also  be  found  for  all  the  bursting  sphere  ^’uns  presented.  This  has  been 
done  and  the  results  plotted  in  Fig.  50.  On  this  graph  the  minimum  sphere 
energy  for  ignition,  is  plotted  as  a function  of  the  non-dimensional 
energy  density,  q. 

It  can  be  seen  that  Fig.  30  bears  a striking  resemblance  to  Fig.  23, 
the  delay  time  plot  for  the  combined  bursting  sphere  runs.  On  Fig.  30, 
the  curves  for  initial  temperature  ratio  equal  to  33,  initial  temperature 
ratio  equal  to  11,  and  initial  temperature  ratio  equal  to  initial  pressure 
ratio  are  drawn.  Moreover,  the  curves  occupy  the  same  relative  positions 
as  those  drawn  on  Fig.  23.  It  is  not  surprising  that  the  two  figures  are 
similar  because  the  delay  time  and  the  bursting  sphere  energy  are  intimately 
related,  as  has  been  discussed  previously.  There  is  one  important  difference 
between  the  two  figures  however.  Fig.  23  shows  that  the  delay  times  for 
the  runs  under  investigation  s]>an  about  two  orders  of  magnitude.  On  the 
other  hand.  Fig.  30  indicates  that  the  minimum  energies,  for  ignition  for 
the  same  runs  span  about  seven  orders  of  magnitude.  Thus  the  minimiun  ignition 
energy  is  a much  strongei  function  of  the  energy  density  than  the  delay  time. 

All  throe  of  the  curves  on  Fig,  .'lO  show  the  expected  behavior,  i.o., 
as  the  energy  density  incx'oasea,  the  minimum  energy  for  ignition  decreases. 
Note  that  as  one  approaches  the  cutoff  energy  density,  which  is  analogous 
to  approaching  the  1000'’K  curve  on  Fig.  15,  the  required  onoi'gy  increases 
extremely  rapidly,  eventually  roucl>ing  an  infinite  value  when  one  encounters 
the  Fig.  IS  cutoff  curve.  On  the  utlw>r  lumd  it  is  seen  that  as  one  approaches 
very  large  energy  densities,  very  little  change  in  the  minimum  bursting 
sphere  onox’gy  is  required  to  proeluce  in  ignition. 

As  in  the  Fig.  23  delay  time  |)lut,  runs  7 and  S are  not  shown  on  the 
lerapuratiu'Q  ratio  equal  to  pressur>.>  ritio  curve,  because  the  incomplete 
numerical  Integration  wfouUl  yield  en\>nemis  data.  Nevertheless,  the  sound 


speed  effects  are  evident  in  this  curve  and  in  the  curve  of  temperature 
ratio  equal  to  11.  The  same  reasoning  discussed  earlier  applies.  That  is, 

■as  the  initial  sphere  temperature  ratio  increases,  so  does  the  sphere 
sound  speed.  This  enables  the  inward  moving  rarefaction  fan  to  release 
the  sphere  energy  to  the  flow  field  more  quickly  and  from  a minimum 
ignition  energy  standpoint  more  efficiently. 

In  summary,  it  should  be  noted  that  a bursting  sphere  may  be  charac- 
terized by  an  instantaneous  addition  of  the  energy  E^.  This  is  because 
at  time  zero  plus  all  of  the  energy  which  is  ever  going  to  be  present  is 
already  present.  We  now  wish  to  investigate  a number  of  flow  fields  in 
which  a finite  amount  of  time  is  required  to  add  the  energy  E^  to  the  high 
pressure  source  region, 

VI.  ENERGY  ADDITION  RESULTS 

To  determine  the  effepts  of  finite  rate  energy  deposition  times  on 

ignition  behavior  in  tho  computer  gomi'ated  flow  field,  a nun)bor  of  computer 

runs  wore  made  with  a non-dimensional  deposition  time  of  t/t^  = .02.  In 

addition  one  run  was  made  with  a nqn-dimenHional  deposition  time  of  t/t^  » .01. 

'IMs  case  will  bo  discussed  last.  The  mechanics  of  how  this  energy  is  added 

14 

with  a finite  rate  in  the  CLOUD  program  is  discussed  in  Adamesyk. 

To  obtain  a cox'tuin  consistency  between  the  finite  rate  onox'gy  addition 
flow  fields  and  those  generated  by  a iiuvsting  sphere,  it  was  decided  to  force 
the  pro.ssurp  and  temperature  x'utios  in  the  source  regions  of  tho  enex'gy 
addition  runs  at  the  instant  the  addition  was  completed  to  be  equal  to  the 
initial  pressure  and  tomporuture  laitios  in  the  source  regions  of  the  cor- 
responding bursting  spheres.  Thus  wi  have  at  the  instant  the  depu.sition  is 
completed  the  same  temperature  ratios  and  iience  the  same  sound  speed  effects 


as  were  encountered  for  the  bursting  sphere  case.  We  must  note  two  things 
however.  First,  since  the  wave  system  begins  propagating  at  non-dimensional 
time  t/t^  = 0+,  by  the  time  the  deposition  is  complete,  the  source  region 
is  not  completely  homogeneous.  Nevertheless  the  inner  portion  of  the 
source  region  is  homogeneous  since  the  inward  travelling  wave  does  not 
reach  this  area  until  after  the  deposition  is  completed.  Second,  the 
addition  is  assumed  to  occur  al  constant  volume.  Again  because  of  the  wave 
system,  this  is  not  strictly  true  but  as  the  change  in  volume  during  this 
time  interval  is  small,  the  effect  is  ignored. 

To  yield  the  proper  conditions  at  t/t^  = .02,  we  must  have  the  proper 
initial  conditions  for  the  energy  addition  runs.  Thi.s  is  determined  as 
follows.  Recall  the  bursting  .sphere  energy  is  given  by, 


_ ••  P{»  \/ 
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In  terms  of  intensive  variables  this  iiay  bo  re-written, 
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whore  w represents  the  mass  in  the  source  region.  Adaraseyk^^  then  defines 
a now  variable  which  Is  the  onorgy  input  term  for  the  finite  rate  deposi- 

tion CLOUD  runs.  is  then  given  by  the  expression, 

Qf  ■-  ■ ■ 
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= tV- 

Q^.  is  then  used  to  characterise  the  energy  input  for  the  CLOUD  progriiw 
where  the  values  of  and  are  those  values  obtained 

for  the  corresponding  bursting  spheri  ca.se. 
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It  remains  to  determine  tl'.e  initial  pressure  and  temperature  ratios 
for  the  energy  addition  cases.  For  all  of  these  runs  the  initial  pressure 
ratio  was  taken  to  be  unity  i.(;.,  ambient  pressure.  However,  the  initial 
temperature  ratio  is  not  necessarily  equal  to  unity.  We  can  solve  for  the 
initial  temperature  ratio  in  the  following  manner.  Since  we  are  assuming 
a constant  volume  energy  addition,  we  may  write: 
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Then  from  the  perfect  gas  law. 
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where  ADDITION  “ '^BURSTING  SI’llKRK  ” *^o  constant 

drops  out  of  the  above  expressions.  Then, 

f^»/T>)iNii(tf<  <AotnT  >0.0  ^ fTa/To^  ivjact>jJ6  imitk' 
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or  since  (f^/PeilNrtfcy  i^lo»T«M  •“■I  froi,i  above,  we  have 
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Then  using  oqns.  37  and  .3d,  eight  energy  addition  runs  wore  made.  Those 
runs  are  displayed  in  Table  8.  For  each  of  these  runs,  the  computer  output 
was  chockod  to  verify  that  the  proper  conditions  were  obtained  at  the  instant 
the  deposition  was  complotud.  la  all  cases,  verification  was  obtained. 
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TABLE  A 


FLOW  FIELD 
RUN 

CORRESPONDING 
BURSTING  SPHERE 
RUN 

PRi-SSURE 

RATIO 

TEMPERATURE 

RATIO 

Qf 

12 

1 

1 

.5 

81.25 

13 

2 

1 

.25 

81.875 

14 

3 

1 

.125 

82.1875 

15 

4 

1 

.0625 

82.34375 

16 

5 

1 

1 

162.5 

17 

6 

I 

1 

327.5 

18 

7 

1 

1 

657.5 

19 

8 

1 

1 

1317.5 

20  * 

1 

1 

• 

81.25 

ic 

Run  20  has  deposition  time,  t/t^  = .01 

All  other  runs  have  deposition  times,  t/t^  =.U2 


The  plots  of  '1'qi^qup  vs.  cell  number,  vs.  cell  number  and 

explosion  loci  are  not  shown  for  the  cnorj}/  addition  runs  as  they  arc 
qualitatively  similar  to  their  bursting  sphere  counterparts.  However, 
it  has  been  worthwhile  to  display  the  minimum  energy  for  ignition  data  for 
the  energy  addition  runs.  Tiiii  has  been  done  in  fig.  51.  On  this  graph 
the  minimum  energy  required  for  ignition  for  u given  energy  addition  run 
has  boon  plotted  v.s,  the  non-dimensional  energy  density,  q of  the  cer- 
rosponding  bursting  sphere  run  Also  shown  on  this  grapli  for  the  sal.e  of 
comparison  are  the  corresponding  bursting  sphere  curves,  taken  from  llg.  50. 
As  before,  data  for  runs  IB  and  10  are  net  shown  because  of  incomplete 
numerical  integrations  for  theie  runs. . 

It  cun  be  seen  on  t'ig.  51  that  fur  most  of  the  range  of  q,  tlie  minimum 
cnorgy  required  to  obtain  an  i-piitioi  in  an  eneigy  addition  flow  field  is 
less  than  the  amount  required  in  a corresponding  bursting  sphere  flow  field. 
A qualitative  e.\pl.inat ion  ot  tii;  bef  ivior  depends  on  two  oliservations. 


first,  in  the  bursting  sphere  cane  il,.‘  initial  compression  wave  is  a .shock. 
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In  the  energy  addition  case,  the  initial  compression  wave  is  an  acoustic 

wave  and  builds  up  to  a shock  as  the  energy  is  being  deposited  in  the 

source  region.  As  it  is  known  that  an  acoustic  wave  will  do  more  compressive 

work  than  a shock  wave,  we  may  conclude  that  the  lead  wave  in  the  energy 

addition  case  will  have  higher  pressures  and  temperatures  in  its  wake  than 

will  the  corresponding  lead  ware  in  the  bursting  sphere  case.  This  is 

conducive  to  larger  values  of  hence  smaller  values  of  minimum 

IS 

ignition  energy  for  the  energy  addition  cases.  Second,  because  ot  the 
time-dependent  nature  of  the  deposition  in  the  energy  addition  cases,  the 
resulting  lead  wave  may  "focus"  certain  particles  of  the  flow  field  to 
produce  local  hot  spots.  This  will  then  also  lead  to  larger  values  of 
’^CLOUU  smaller  values  of  mi.uimum  ignition  energy. 

At  the  far  left  of  Fig.  51,  where  the  values  of  q are  small,  the  energy 
addition  curves  arc  seen  to  merge  int('  the  bursting  sphere  curvo.s.  What 
is  happening  here  i.s  that  at  tl»ese  low  energy  densities,  ignition  boeoiuos 
hard  to  produce  as  we  are  approaching  the  eriucul  cutoff  curve  of  Fig.  15. 
Thus  both  the  energy  addition  and  bursting  sphere  curves  approach  a vortical 
asymptote  and  because  of  this  are  sti>n  tu  merge  into  each  other.  It 
sliould  also  be  noted  that  at  q - .TtiU  rl\e  bursting  sphere  curve,  for  the 
temperature  ratio  equal  to  .>5  is  seen  tu  intersect  Us  corresponding  energy 
addition  curve,  it  is  l>elieved  that  ilris  Is  cn.sevi  by  an  inceiaplote 
numerical  integration  In  run  U and  tlms  not  ro|resentativc  of  true  physical 
behavior.  However,  this  has  not  been  cuntirmed  at  tiu‘  time  of  this  writing. 

Finally,  we  may  e.\amine  one  further  piienomv  urn.  Wo  can  define  a new 
variable  p,  called  llie  juiwer  d.'nsity  »>  the  following  relation. 


r>fcpo^tttON 


where  q has  been  defined  previoasly  aiul  represents 

the  non-dimensional  energy  deposition  cime.  We  may  look  at  bursting  sphere 
run  1 and  energy  addition  runs  12  and  10.  1-or  each  of  these  case,  q = 65. 
Run  12  has  a deposition  time  of  ==  .02  whereas  run  20  has  a deposition 
time  of  t/t^  = .01,  Because  run  1 is  a bursting  sphere  it  is  known  to  have 
a deposition  time  t/t^  = 0.  For  these  cases  then  the  power  densities  have 

’•i 

boon  determined  from  eqn.  39.  The  inverses  of  these  power  densities  have 
also  been  determined.  These  data  are  shown  in  Table  9. 
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Fig.  52  is  u plot  of  the  iiiiumum  energy  required  i'or  ignition  vs.  the 
quantity  l/p  for  the  above  three  c.ase>’.  ibo  data  points  are  connected  by 
stx'aight  dashed  lines  to  inUieuto  tluii  the  exact  shape  of  the  curve  i.s  not 
known  at  this  time.  However,  "'rum  tlie  three  data  points  plotted,  we  can 
say  that  the  curve  appears  to  '.'xhibit  relative  w.\treinum  behavior.  This 
.suggests  that  there  may  possibly  be  optimal  power  densitie.s  for  ignition. 
This  has  not  been  investigated  furihe-  at  the  time  of  writing. 


VU  C0Ni:i.USl0.N.S 


We  have  .seen  that  tlu*  ignition  li'ltavinr  in  a non-ste.uly,  mm-ldeui  flovs' 
field  is  a complicated  function  of  th.‘  flow  variables,  particularly  cell 
energy  or  temperature  and  cell  speci'tc  volume.  These  quantities  are 
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themselves  complicated  functions  of  the  initial  conditions  in  the  source 
region. 

We  have  also  seen  that  the  important  quantities  to  determine  are  the 
delay  times  to  ignition  and  the  minimum  energies  required  for  ignition,  and 
that  both  of  these  quantities  tend  to  decrease  as  the  non-dimensional 
energy  density  increases. 

We  have  also  studied  the  differences  in  ignition  behavior  as  the 
source  was  changed  from  a bursting  sphere  to  a finite-rate  energy  deposi- 
tion source,  and  found  that  for  most  values  of  i|,  the  finite  rate  energy 
addition  cases  required  a smaller  amount  of  dimensional  energy  to 
produce  an  ignition  in  the  flov>  field. 

Finally  we  have  briefly. oamincd  the  effects  of  power  density  on 
minimum  energy  required  for  ignition  and  sl\ow  that  there  may  be  optimal 
power  densities  to  pi'oduce  an  ignition, 

VI H.  SUGGESTIONS  FOR  FUTURE  RliSEARCil 

A number  of  variables  have  been  hold  constant  during  this  investigation. 
Tho.so  could  bo  varied  to  yield  additional  data.  For  c.\umple  one  could  change 
the  temperature  cutoff  limit  from  10Ul>®K  to  12U0  “K  to  see  what  effect 
it  has  on  tho  ignition  bohavior  in  the  resulting  flow  fields.  One  could 
also  substitute  the  values  for  tho  ac  ivuiion  energy  and  pro-exponential 
factor  from  tho  Schott  and  Kiiney  data  for  tho  \uluos  obtained  from  Strohlow. 
One  could  also  use  source  region.^;  which  contain  for  e.vample  helium 
(Yj  = 1.667).  In  addition  it  would  V"  Instructive  to  obtain  more  data 
points  for  Fig.  32  at  different  power  densities. 

There  are  tiu’ue  limitations  of  the  model  used  in  this  Investigation 
which  might  he  removed  in  future  worh.  first  , ;i  more  detailed  reaction 


kinetics  scheme  could  be  incorporated.  Second,  when  a particle  ignites, 
it  then  adds  heat  to  the  neighboring  particles,  thus  altering  the  flow 
field.  This  heat  addition  could  also  be  included  in  the  model.  Third, 
the  combustible  mixture,  which  was  assumed  to  have  the  properties  of  air 
in  this  investigation  could  be  replaced  by  a more  realistic  combustible 
mixture. . 

Finally,  the  point  source  solution  for  spherical  blast  waves  could 
be  used  as  a means  of  establishing  a flow  field  for  ignition  behavior 
studies.  As  the  point  source  solution  provides  an  ideal  flow  field  it 
would  contrast  nicely  with  the  non-ideal  bursting  sphere  and  energy 
addition  flow  fields  studied  in  this  report. 
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